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Abstract 

Inclusive neutrino-nucleus cross sections are calculated using a consistent relativistic mean-field 
theoretical framework. The weak lepton-hadron interaction is expressed in the standard current- 
current form, the nuclear ground state is described with the relativistic Hartree-Bogoliubov model, 
and the relevant transitions to excited nuclear states are calculated in the relativistic quasiparti- 
cle random phase approximation. Illustrative test calculations are performed for charged-current 
neutrino reactions on ^^C, ^^O, ^^Fe, and ^'^'^Pb, and results compared with previous studies and 
available data. Using the experimental neutrino fluxes, the averaged cross sections are evaluated 
for nuclei of interest for neutrino detectors. We analyze the total neutrino-nucleus cross sections, 
and the evolution of the contribution of the different multipole excitations as a function of neutrino 
energy. The cross sections for reactions of supernova neutrinos on ^^O and ^^^Pb target nuclei are 
analyzed as functions of the temperature and chemical potential. 

PACS numbers: 21.30.Fe, 21.60.Jz, 23.40.Bw, 25.30.-c 
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I. INTRODUCTION 



Neutrino-nucleus reactions at low energies play an important role in many phenomena 
in nuclear and particle physics, as well as astrophysics. These reactions present extremely 
subtle physical processes, not only because they involve the weak interaction, but also be- 
cause they are very sensitive to the structure of nuclear ground states and excitations, 
i.e. to the solution of the nuclear many-body problem that includes the strong and elec- 
tromagnetic interactions. The use of microscopic nuclear structure models in a consistent 
theoretical framework is therefore essential for a quantitative description of neutrino-nucleus 
reactions [l^ . Detailed predictions of neutrino-nucleus cross sections are crucial for the inter- 
pretation of neutrino experiments, detection of neutrinos produced in supernova explosions 
and understanding the underlying nature of these explosions Neutrino-nucleus reactions 
which occur in a type II supernova could also contribute to the nucleosynthesis js, 4], but 
more data on cross sections are necessary for a more complete understanding of this process, 
as well as the supernova dynamics. 

n n 

Data on neutrino-nucleus cross sections have been _obtained by the LSND [5|, 16[ and 



10 



Ul, but only for and ^^Fe 



KARMEN 0, y, 3 collaborations, and at LAMPF 
target nuclei. New experimental programs are being planned which will provide essential 
data on the neutrino-nucleus reactions, and also help to improve the reliability of present 
cross-section calculations. These include the spallation neutron source (SNS) at ORNL, 
where neutrinos produced by pion decay at rest will enable measurements of cross sections 
for a wide range of target nuclei 12, ll3| , and the promising "beta-beams" method for the 
production of pure electron neutrino-beams by using the /5-decay of boosted radioactive 
ions Q, [l^. Cross sections for neutrino energies in the range of tens of MeV could be 
measured, and these reactions are particularly interesting for supernova studies [l^. 



Weak interaction rates at low energies have been analyzed employing a 



scopic approaches, principally in the 



phase approximation (RPA) [20 



21 



rameworks of the shell model 
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22 



hybrid models of CRPA and shell model 
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variety of micro- 
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9L the random 
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3I32 



33|. The 



continuum RPA (CRPA) ^ 
29|, |30| , and the Fermi gas model 
shell model provides a very accurate description of ground state wave functions. The descrip- 
tion of high-lying excitations, however, necessitates the use of large model spaces and this 



often leads to computational difficulties, making the approach applicable essentially only to 
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light and medium-mass nuclei. For systematic studies of weak interaction rates throughout 
the nuclide chart, microscopic calculations must therefore be performed using models based 
on the RPA. 

Hybrid models combine the shell-model and CRPA in such a way that occupation prob- 
abilities of single-particle states in a specific nucleus are determined by shell model cal- 



culations, and then inserted into the CRPA 



2^. In general the CRPA employs different 



interactions for the calculation of the nuclear ground state (for instance, the Woods-Saxon 
potential), and in the residual CRPA interaction (G matrix from the Bonn potential, or 
the Landau-Migdal interaction), and thus additional parameters are required in order to 
adjust the calculated rates to data. A more consistent approach, based on the quasipar- 
ticle RPA with Skyrme effective interactions, has been employed in calculations of weak 



interaction rates 



221 ]. Although in this framework the residual RPA interaction is derived 



from the same energy functional which determines the nuclear ground state, some terms are 
usually omitted, and pairing correlations require the adjustment of additional factors. A 
fully consistent theoretical framework for the description of charge-exchange excitations in 
open shell nuclei, based on Skyrme effective interactions, has only recently been developed 



3J], but not yet employed in the analysis of neutrino- nucleus reactions. Neutrino-nucleus 



cross sections have also been calculated using the ab-initio no-core shell model based on 



realistic NN and three-body interactions 



35| , and the shell model with an improved Hamil- 



tonian which properly takes into account the spin-isospin interactions 36|] . The importance 
of improved calculations of neutrino-nucleus cross sections for neutrino-oscillation studies 

n 

has been demonstrated in the recent reanalysis of the LSND experiment using the 
particle- number projected quasiparticle RPA [38i], which has shown an enhancement of the 
neutrino-oscillation probability when compared to previous studies. 

Although the general expressions for the transition matrix elements relevant for the calcu- 
lation of neutrino-nucleus cross sections have been known since many years 3^, I40], it is only 
more recently that systematic calculations have been performed in open-shell nuclei by mak- 
ing use of modern effective interactions in the description of both nuclear ground states and 



excitations 



221 ]. Reliable prediction of weak interaction rates in nuclei necessitates a fully 



consistent description of the structure of ground states and multipole excitations. Among 
the relevant charge-exchange excitations, the isobaric analog state (IAS) and Gamow- Teller 
resonance (GTR) have been the subject of extensive experimental and theoretical studies. 
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Much more limited are the data and theoretical predictions for properties of excitations of 
higher multipolarities at finite momentum transfer. 

In this work we analyze charged-current neutrino-nucleus reactions by employing a fully 
consistent microscopic approach based on relativistic energy density functionals, and also 
including pairing correlations in the description of open-shell target nuclei. An essential 
advantage over most current approaches is the use of a single universal effective interaction 
in calculations of both ground-state properties and multipole excitations of nuclei in various 
mass regions of the chart of nuclides. 

Of particular interest for the present study are rates for neutrino-nucleus reactions in 
the low-energy range below 100 MeV, which play an important role in many astrophysical 
processes, including stellar nucleosynthesis. A quantitative description of nucleosynthesis 
of heavy elements during the r-process necessitates accurate predictions of neutrino-nucleus 
cross sections not only in stable nuclei, but also in nuclei away from the valley of /^-stability. 
Since nuclei are used as detectors for solar and supernovae neutrinos, as well as in neutrino 
oscillation experiments, it is important to describe the neutrino detector response in a con- 
sistent and fully microscopic theory. Finally, a quantitative estimate of neutrino-nucleus 
reaction rates will provide information relevant for feasibility studies and simulations of a 
low-energy beta beam facility, which could be used to produce neutrino beams of interest 
for particle physics, nuclear physics and astrophysics jl5| . 

In Sec. [Til we outline the basic formalism used in the evaluation of neutrino- nucleus cross 
sections. In Sec. IIIII the relativistic Hartree-Bogoliubov model (RHB), and the proton- 
neutron relativistic quasiparticle random phase approximation (PN-RQRPA) are briefly re- 
viewed. Section [IV] includes several test cases, and the calculated cross sections are compared 
with results of previous theoretical studies and the available data. Cross sections for su- 
pernova neutrinos are analyzed in Sec. |Vl and finally Sec. IVII summarizes the results of the 
present investigation and ends with an outlook for future studies. 

II. NEUTRINO-NUCLEUS CROSS SECTIONS 

In the present study we consider the charged-current neutrino-nucleus reactions: 

t^i +z Xn — ^N-l + ^ ) (1) 
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where / denotes the charged lepton (electron, muon). Detailed expressions for the reaction 



rates and the transition matrix elements can be found in Refs. 
neutrino-nucleus cross section reads 



31 



391]. The charged-current 



spins 

where pi and Ei are the momentum and energy of the outgoing lepton, respectively. The 
Hamiltonian Hw of the weak interaction is expressed in the standard current-current form, 
i.e. in terms of the nucleon J\{x) and lepton i\{x) currents 

Hw = / dxJx{x)3\x) , (3) 

and the transition matrix elements read 

G_ 

The multipole expansion of the leptonic matrix e 
structure for the nuclear transition matrix elements 
nucleus cross section. In the extreme relativistic limit — (ERL), in which the energy of the 
outgoing lepton is considered much larger than its rest mass, the differential neutrino-nucleus 
cross section takes the form 

/da A _ 2Glcos^ec Ef 



{f\Hw\i) = -^Ix J dxe'''''^{f\j\x)\z) . (4) 

ement Zac"**^^ determines the operator 



3l| , and the expression for the neutrino- 



2q2 2 2 



j>i 



z \l q / z 
^ J>0 1^1 ) 



(5) 



where is the Fermi constant for the weak interaction, 6c is the Cabbibo's angle, 6 
denotes the angle between the incoming and outgoing leptons, the energy of the lepton 
in the final state is Ei, and the 4- momentum transfer is g = {qo, q). The nuclear transition 
matrix elements between the initial state | Jj) and final state \Jf), correspond to the charge 
Aij, longitudinal Cj, transverse electric Tf^, and transverse magnetic T^^^'^ multipole 



5 



operators. These are expressed in terms of spherical Bessel functions, spherical harmonics, 
and vector spherical harmonics: 

Mf ^jj{Kx)Yf{n,) (6) 

Mfi^jL{^x)Yfl,{Q^), (7) 
where we use the standard notation n = \q\. There are four transition operators in Eq. ([5]): 



the Coulomb multipole operator 



MjMix) = F,^Mf{x)~t 



rriN 



-(Fa - qoFp)J:r{x 



(8) 



the longitudinal operator 
^jm{x) = 

the transverse electric operator 



' 1 ^ ■ 
rA — 

2 rriAT 



^Tix) , 



"^flvAx) 



K 



rriN 



+ iFA^'f{x) 



(9) 



(10) 



and the transverse magnetic operator 



K 



rriN 



+ FA^fix) 



(11) 



where all the form factors are functions of q^. These multipole operators contain seven 
fundamental operators expressed in terms of spherical Bessel functions, spherical harmonics, 
and vector spherical harmonics 3^: 



Afix) = MrAx)-V 



Mfix) 



rM , 



K 



S"f fa;) 



J 



2J + r 



M 



M 

JJ+1 



X 



I J+l 

2J+1 



M'jUix) 



J:y{x) = MjAx)a 



(12) 
(13) 

(14) 
(15) 
(16) 

(17) 
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VL^{x) = Mf{x)cT-V. (18) 

By assuming conserved vector current (CVC), the standard set of form factors reads 331. l39| : 

i^iV)=fl+ ... V (19) 



(855Me\/)2 



/ ^2 



-2 



/ 2 \ -2 

^ _ . (22) 

The neutrino- nucleus cross section is evaluated using Eq. (j5]), with the transition matrix 
elements between the initial and final states determined in a fully microscopic theoretical 
framework based on the Relativistic Hartree-Bogoliubov (RHB) model for the nuclear ground 
state, and excited states are calculated using the relativistic quasiparticle random phase 
approximation (RQRPA). More details about the RHB+RQRPA framework are presented 
in the next section. 



III. THE PROTON-NEUTRON RELATIVISTIC QUASIPARTICLE RANDOM 
PHASE APPROXIMATION 

Here we outline the essential features of the theoretical framework that will be used to 
evaluate inclusive neutrino-nucleus cross sections: the proton-neutron (PN)-RQRPA based 
on relativistic Hartree-Bogoliubov model. The detailed formalism, as well as recent appli- 
cations of these models in the description of exotic structure phenomena in nuclei far from 



stability, have been reviewed in Refs. 



41 



42|. 



The RQRPA has been formulated in the canonical single-nucleon basis of the relativistic 
Hartree-Bogoliubov (RHB) model in Ref. {43}], and extended to the description of charge- 
exchange excitations (proton-neutron RQRPA) in Ref. [4^. In the RHB framework the 
ground state of a nucleus can be written either in the quasiparticle basis as a product of 
independent quasi-particle states, or in the canonical basis as a highly correlated BCS-state. 
By definition, the canonical basis diagonalizes the density matrix and is always localized. 
It describes both the bound states and the positive-energy single-particle continuum. The 



QRPA equations in the canonical basis include the matrix elements V^v^'a of the residual 
particle-hole {ph) interaction, and the matrix elements V^^'aa' of the particle-particle (pp) 
(pairing) interaction, as well as certain combinations of occupation factors u^, of the 
canonical single-nucleon states. In addition to the configurations built from two-quasiparticle 
states of positive energy, the relativistic QRPA configuration space must also include pair- 
configurations formed from the fully or partially occupied states of positive energy and the 
empty negative-energy states from the Dirac sea. The inclusion of configurations built from 
occupied positive-energy states and empty negative-energy states is essential for current 
conservation, decoupling of spurious states, and a quantitative description of excitation 



energies of giant resonances 



43|. 



The RHB+RQRPA model is fully consistent: in the particle- hole channel effective La- 
grangians with density-dependent meson-nucleon couplings are employed, and pairing corre- 



45j. Both in 



lations are described by the pairing part of the finite range Gogny interaction 
the ph and pp channels, the same interactions are used in the RHB equations that determine 
the canonical quasiparticle basis, and in the matrix equations of the RQRPA. This is very 
important because the energy weighted sum rules are fulfilled only if the pairing interaction 
is consistently included both in the static RHB and in the dynamical RQRPA calculation. 
In the present work all RHB plus QRPA calculations have been performed using one of the 
modern and most accurate meson-exchange density-dependent relativistic mean-field effec- 
tive interactions - DD-ME2 [s^ in the ph channel, and the finite range Gogny interaction 
DIS [45^ in the pp channel. 

Charge-exchange excitations have been analyzed with the proton-neutron RQRPA using 



effective Lagrangians characterized by density-dependent meson-nucleon couplings 



4^. The 



RHB plus PN- RQRPA framework has been successfully employed in studies of isobaric ana- 



log resonances, Gamow- Teller resonances 



44| . the relation between the excitation energies 



of IAS and GTR and the evolution of neutron skin in neutron- rich nuclei |46|, and /3-decay 
rates of r-process nuclei 47|, |48|. The model includes both the T = 1 and T = pair- 
ing channels, and represents a relativistic extension of the fully-consistent proton-neutron 
QRPA introduced in Ref. [49 1. 
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The matrix equations of the PN-RQRPA read 



^ A B 




E 



X 



1 
-1 



B A 

with the matrices A and B defined in the canonical basis 



(23) 



+ {UpVnUpiVn' + VpUnVp'U + {UpUnUpiUn' + VpVnVp'Vn') KfZ,,„, 



Bpn,p'n' = {UpVnVpiUn' + VpUnUp'Vn')V^p,^^, 

- iUpUnVp'Vn' + VpVnUp'Un') Vp^p,„, . (24) 

The proton and neutron quasiparticle canonical states are denoted by p, p', and n, n', 
respectively. V^^ is the proton-neutron particle-hole residual interaction, and V^^ denotes 
the particle-particle interaction. The canonical basis diagonalizes the density matrix, and 
the corresponding eigenvalues are the occupation probabilities fp(„). Since on the other hand 
the canonical basis does not diagonalize the Dirac single- nucleon mean-field Hamiltonian ho, 
nor the pairing field A, the off-diagonal matrix elements H^\, and Hpp, appear in Eq. (Ell): 

Hll, = {u^u^i - v^v^i)h^^, - {u^v^' + v^u^>)A^^> . (25) 

The solution of the RQRPA matrix equation (!23ll determines the excitation energies Ex and 
the corresponding forward- and backward-going amplitudes and Y''^, respectively. In the 
evaluation of the neutrino- nucleus cross section (Eq. ([5])), for each transition operator O"^ 
the matrix elements between the ground state of the even-even (N,Z) target nucleus and the 
final state are expressed in terms of single-particle matrix elements between quasiparticle 
canonical states, occupation probabilities, and the corresponding RQRPA amplitudes: 

{JfWO'm = EiPWO'W^) (xinUpVn - YpiVpUr,) . (26) 
pn 

In the calculations performed in this work we have considered transitions between the lO"*") 
ground state of a spherical even-even target nucleus and excited states in the odd-odd nucleus 
with multipolarities JJ = 0^ - 7^. 

The spin-isospin-dependent interaction terms are generated by the exchange of it and p 
mesons. Although the direct one-pion contribution to the nuclear ground state vanishes at 
the mean-field level because of parity conservation, the pion nevertheless must be included 



in the calculation of spin-isospin excitations which contribute to the neutrino-nucleus cross 
section. The corresponding particle-hole residual interaction of the PN-RQRPA is derived 
from the Lagrangian density 

C:i^ = -g^^pYp^T^ - A^^^Yd.nri^ , (27) 

777, 7r 

where arrows denote vectors in isospin space. The residual two-body interaction reads 

- (—] rir2(SiVi)(5]2V2)/^.(ri,r2) , (28) 
where -Dp(^) is the meson propagator 

1 g-™p(ir)kl-r2 

^pW = 1 r- ' (29) 

^ ' 47r |ri — r2\ 

and 

/ rr n \ 

(30) 

For the p-meson density-dependent coupling strength we choose the functional form used in 
the DD-ME2 effective interaction 

9p{Pv) = 9p{Psat)e-''^^''-'^ , (31) 

with X = Pv/Psat, and Psat is the nucleon density at saturation in symmetric nuclear matter. 

For the pseudovector pion-nucleon coupling we use the standard values: /^/47r = 0.08, 
and = 138 MeV. The derivative type of the pion-nucleon coupling necessitates the 
inclusion of the zero-range Landau-Migdal term, which accounts for the contact part of the 
nucleon-nucleon interaction 

V5. = g'(—] riTaSi ■ S25(ri - r2) , (32) 

with the parameter g' adjusted in such a way that the PN-RQRPA reproduces the exper- 




imental values of GTR excitation energies 4J|. For the DD-ME2 effective interaction the 
value of g'=0.52 has been adjusted to the position of the GTR in ^°^Pb. 



In most applications of the RHB and RQRPA models 
Gogny force has been used in the the pp-channel: 



41 



421 ]. the pairing part of the 



\/pp(l,2) = ^ e-[(''i-'"2)/'''l' (IVi + BiP'' - HiP^ - MiP"P^), (33) 

=1,2 
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with the set DIS [45|] for the parameters /Xj, Wi, Bi, Hi and Mj {i = 1,2). This force has 
been very carefully adjusted to pairing properties of finite nuclei all over the periodic table. 
In particular, the basic advantage of the Gogny force is the finite range, which automatically 
guarantees a proper cut-off in momentum space. 

IV. ILLUSTRATIVE CALCULATIONS: CROSS SECTIONS FOR NEUTRINO 
DETECTOR RESPONSE 

The theoretical framework described in the previous two sections can systematically be 
applied in studies of charged- current neutrino reaction rates with target nuclei of arbitrary 
mass. We have performed several illustrative calculations of cross sections for reactions with 
nuclei of interest for neutrino detector response: ^^C, ^^O, ^^Fe, and ^^'^Pb. The inclusive 
cross sections, which sum the contributions from transitions to all possible final states, are 
given as functions of neutrino energy, and can be averaged over the experimental neutrino 



flux when available, e.g. from the decay at rest of yU"*", and the decay in flight of vr 



aa 
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111 ]. The fiux-averaged cross sections provide a crucial test for the validity of the 
theoretical approach used for modeling neutrino-nucleus reactions. 

For charged- current reactions the cross section Eq. ([2]) must be corrected for the distortion 
of the outgoing lepton wave function by the Coulomb field of the daughter nucleus. In order 
to be able to compare our results with previous studies ^^e the same prescription 

for the Coulomb correction. The cross section can either be multiplied by a Fermi function 
obtained from the numerical solution of the Dirac equation for an extended nuclear charge 

:t of the Coulomb field can be described by 



51 



52 



531 ]. in which the lepton momentum 



distribution [SOj] or, at higher energies, the effect 
the effective momentum approximation (EMA) 
pi and energy Ei are modified by a constant electrostatic potential evaluated at the center 
of the nucleus: 



pff = ^jEt^f - (34) 

Ef// = Ei- VciO), (35) 

where Vc{0) = —3Za/{2rc), the nuclear charge radius is denoted by rc, and m/ is the mass of 
the outgoing lepton. In calculations with EMA the original lepton momentum pi and energy 
El appearing in the expression for the cross section, are replaced by the above effective 
quantities. 
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The two prescriptions for the Coulomb correction are tested in the calculation of the 
inclusive cross section for neutrino scattering on a high-Z nucleus: 2°'^Pb(z/e, e~)2°^Bi. In 
Fig. [T] we display the corresponding cross sections corrected either with the Fermi function, 
or calculated with the EMA. As already shown in previous studies, the Fermi function 
correction overestimates the cross sections at higher neutrino energies, where the EMA 
provides a more reliable correction. Therefore, following the prescription from Ref. {22], we 
will use the Fermi function correction in the range of neutrino energies for which the cross 
section is below the corresponding EMA value, whereas the EMA will be employed at higher 
energies (cf. the solid curve in Fig. [1]). 

One of the most extensively studied neutrino-nucleus reactions is ^2C(z/e, e~)^2N. This 
reaction is particularly important because ^^C is used in liquid scintillator detectors, and data 
on the cross section are available from the LSND 0,3, and KARMENQ,y, collaborations, 
as well as from LAMPF 10|,lll(|- In order to illustrate the contributions of different multipole 
excitations, in Fig. [2] we plot the inclusive cross section for the '^'^C{i'e,e~y'^N reaction as 
function of the neutrino energy, obtained by successively increasing the maximal allowed 
angular momentum in the sum over J in Eq. ([5]): from Jmax = 0^ to Jmax = 7^- One notices 
that the largest contributions arise from J = 1"^ and J = 2^, and that the contribution of 
higher multipolarities gradually decreases. In fact, in this figure one can not distinguish the 
cross sections calculated with Jmax = 6^ and Jmax = 7^? for the whole interval of neutrino 
energies. 

The results of model calculations can be compared with available data by averaging the 
cross section over the neutrino flux f{Ei,), which depends on the specific neutrino source 

, , !dE„a,(E,)f(E,) 



For Ve the Michel flux from the decay at rest (DAR) of /i"*" is used lOj 



/(i^.e) = ^ K - 2i5;.J , (37) 



whereas for we employ the polynomial fit to the experimental flux obtained from the decay 
in flight (DIE) of 7r+ js^]- In Table [T] the values of the muon and electron neutrino-nucleus 
cross sections for the reactions ^'^C{veiG~Y'^^ and ^"^C^u^^, calculated with the PN- 

RQRPA and averaged over the empirical neutrino fluxes, are compared with the results 
of previous theoretical studies based on the non-relativistic random phase approximation 
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and the shel 



LSND 0, 154, 



model, and with the experimental values: LAMPF lOj, KARMEN 



and 



551]. For the {v^i^i ) reaction the PN-RQRPA cross section is in excellent 



2^, and the CRPA value of Ref. |29|, but 



agreement with the (Q)RPA results of Ref. 
somewhat larger than the experimental cross section. The PN-RQRPA cross section for 
the ^^C(z/e, e^)^^N reaction is in bette r ag reement with experiment when compared with the 
results of the QRPA calculations 0, 38] . We notice that a very good agreement with data 



56j, and with the RPA 



has been obtained in the shell-model calculations of Refs. 
based on the Skyrme interactions SGII and SIII [2^. 

In Fig. [3] the calculated inclusive cross sections are plotted up to 100 MeV neutrino 
energy for three reactions including also heavier nuclei: ^^0(z/e, e~)^^F, ^^Fe(i/e, e~)^^Co, 
and ^°^Pb(z/e, e~)^*'^Bi. The results of the RHB plus PN-RQRPA calculations are compared 
with those of a very recent analysis performed with the nonrelativistic QRPA and using the 
Skyrme interaction SIII 57|]. First we notice the pronounced enhancement of cross sections 
for neutrino reactions on heavier nuclei: from ^^O to ^°^Pb the reaction cross sections increase 
by more than two orders of magnitude. For the reaction on ^^O the present results are in 
excellent agreement with the QRPA - SIII cross sections 57|] at all neutrino energies. This is 
also the case for the reactions on ^^Fe and ^osp]^ ^t neutrino energies above 40 MeV, whereas 
at lower energies the RHB plus RQRPA cross sections appear somewhat below the QRPA 
results. 

In table [III we compare various theoretical results of flux-averaged cross sections for the 
electron neutrino reactions with ^^O, ^^Fe, and ^osp^ target nuclei. For the neutrino flux the 



DAR spectrum Eq. ( |3711 is used, and values obtained with the shell model j56|], the tensor 
model 59], Tamm-Dancoff [60], and various RPA-based calculations are included. For the 
^^O target, the present PN-RQRPA (DD-ME2) cross section agrees with the shell model and 
(Q)RPA results. For heavier targets the differences between various calculations are more 
pronounced. Particularly important is the reaction ^^Fe(z/e, e~)^^Co, because it is the only 
case in which data are available for a medium- heavy target nucleus ol, IgiI. The two self- 
consistent calculations: the QRPA [52] and the present RHB plus PN-RQRPA, which do not 
require any additional adjustments of the model parameters to the specific target nucleus, 
predict values for the cross section: 352 xlO~'^^cm^ (QRPA [57]) and 140 x 10"'^^cm^ (PN- 
RQRPA) which differ considerably, but are still within the uncertainties of the experimental 
result (256±108±43)xl0-^2cjn2_ 
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Obviously various theoretical approaches differ significantly in the predicted neutrino- 
nucleus cross sections, and this will require more detailed studies of the underlying nuclear 
structure that contributes to the neutrino reaction rates. On the other hand, the only 
data for a medium-heavy target nucleus are from the KARMEN collaboration 9|, and this 
result has not yet been confirmed by independent measurement. One hopes that future 
experiments will provide additional constraints for theoretical models that are used in the 
description of weak interaction rates. At the same time the nuclear structure input for 
neutrino-nucleus cross sections must be further analyzed, and checked by improving the 
description of nuclear ground-state properties and spin-isospin excitations relevant for the 
calculation of the reaction rates. 

Next we consider in more detail the contributions from excitations of various multipo- 
larities to the neutrino-nucleus cross sections evaluated within the RHB plus PN-RQRPA 
framework. Fig. H] displays the calculated cross sections for the reaction ^^Fe(z/e, e~)^^Co. 
The solid curve represents the inclusive cross section as a function of the neutrino energy 
calculated by including all multipole transitions with J < 7^, whereas the dashed curve 
corresponds to the cross section evaluated only for transitions to the isobaric analog state 
(IAS) and Gamow- Teller (GT) resonance. The difference between these two cases is repre- 
sented by the dotted curve which, as one notes, rises sharply with energy. Below an incident 
neutrino energy of E^^ ^ 30 MeV the dominant contribution to the cross section originates 
from the two basic charge exchange modes: IAS and GTR. With increasing E^^, however, 
contributions from other multipole transitions become even more important. This resu 



clearly shows that models which take into account only the IAS and GTR transitions 62| . 



cannot provide quantitative predictions of cross sections at neutrino energies E^,^^ > 30 MeV. 

The contribution of different multipole transitions to the cross section for the reaction 
^^Fe(z/e, e^)^^Co, at increasing neutrino energies, is illustrated in Fig. |5l For the neutrino 
energies E^^ = 20, 40, 60, and 80 MeV, we display the partial contributions to the total cross 
section of multipoles from J = 0^ to J = 5"^. At E^^ = 20 MeV the reaction is dominated 
by transitions to the IAS and GTR. At higher energies, however, other multipolarities start 
to play an important role and one notices, in particular, the dominant contribution of the 1~ 
transitions at -Ej,^=60 and 80 MeV. This distribution of multipole transitions, calculated in 
the RHB plus PN-RQRPA model, is similar to the one obtained in the QRPA with Skyrme 
effective interactions 57(1 . The differences in the contributions of higher multipoles can be 
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attributed to the different effective interactions employed in the two models. 



V. CROSS SECTIONS FOR SUPERNOVA NEUTRINOS 

An important application of microscopic models of neutrino-nucleus reactions is the cal- 
culation of cross sections for supernova neutrinos. Accurate modeling of reaction rates on 
nuclei that can be used as targets for the supernova neutrino detectors is, of course, essential 
for studies of supernova dynamics and, in particular, of weak interaction processes which 
determine the evolution of a supernova explosion. In this section we compare the results of 
illustrative RHB plus PN-RQRPA calculations of charged- current reaction rates with those 
obtained in the shell model 17|, and with the CRPA 30(]. The supernova neutrino flux is 
usually described by the Fermi-Dirac spectrum: 

1 

^^^"^ = f^exp[iE,/T]-a] + l ' ^^^^ 

At typical supernova neutrino energies one expects the total cross section for the charged 
current reaction (i/e,e~) to be dominated by the allowed transitions to the IAS and the 
Gamow- Teller resonance states in the daughter nucleus. In Fig. [HI we display the flux- 
averaged cross sections for the reaction ^^0(z/e, e~)^^F, evaluated at different temperatures 
in the interval T = 2 — 10 MeV, and for the chemical potential a=0. The results are in a 
remarkable agreement with those obtained in the shell-model calculation of Ref . 17|] . This 
is a very interesting result, because the two models are based on entirely different concepts, 
and use different effective interactions. An important advantage of the present approach, 
however, is that it can easily be extended to heavier nuclei and to systems far from stability, 
and therefore also applied in the description of the r-process nucleosynthesis. 

We have further tested the model in the case of the heavy target: ^°^Pb(i^e, e~)^°^Bi. For 



the same neutrino flux as in the previous example, the corresponding cross sections are s^ 



nown 



in Fig. [7] and compared, at T = 6, 8 and 10 MeV, with the results from Refs. 2a, 13Q|, l63 |. 
Since at higher temperatures the neutrino flux extends toward higher neutrino energies, the 
calculated cross sections display a pronounced enhancement with temperature (cf. Sec. IIVI) . 
In comparison with the RPA-based results from Refs. 0, present PN-RQRPA 

calculation predicts a slightly lower value for the fiux-averaged cross section at T = 6 MeV, 
whereas an excellent agreement is obtained at T = 8 and 10 MeV. Note that the calculation 
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of Ref . 



631], which includes the IAS, the GT, and the first-forbidden transitions computed in 



the Goldhaber- Teller model, consistently predicts higher values for the flux-averaged cross 
section at all considered temperatures. As it has already been discussed in Ref. [23|, a 
possible origin of this difference could be the use of the Fermi function for the Coulomb 
correction at higher neutrino energies, instead of the EMA. 

Although most calculations of neutrino-nucleus cross sections consider only the case when 
the chemical potential a=0, in astrophysical applications it might be important to perform 
studies of reaction rates for other values of a 30|. We have therefore explored the sen- 



sitivity of the neutrino reaction cross sections to a. In Fig. [8] the contour plot is shown 
for the ^^0(z/e) e~)^^F cross sections as functions of T and a in the Fermi-Dirac spectrum of 
Eq. (1551) . The curves in the plot correspond to constant values of the cross sections. At lower 
temperatures the cross sections depend only very weakly on a. However, already above T=2 
MeV, the cross sections increase much faster for higher values of the chemical potential a. 



VI. CONCLUDING REMARKS 



Detailed microscopic calculations of charged-current and neutral-current neutrino-nucleus 
reaction rates are of crucial importance for models of neutrino oscillations, detection of su- 
pernova neutrinos, and studies of the r-process nucleosynthesis. In this work we have in- 
troduced a consistent microscopic method for calculating neutrino-nucleus cross sections, 
based on the relativistic mean-field framework. In this approach the ground state of a nu- 
cleus is described with the relativistic Hartree-Bogoliubov model, and the neutrino-induced 
transitions to excited nuclear states are computed in the relativistic quasiparticle random 
phase approximation. Since it is based on the self-consistent mean-field approach to nuclear 
structure, the model can be applied to neutrino reactions with target nuclei of arbitrary 
mass throughout the chart of nuclides. By employing universal effective interactions, with 
parameters adjusted to global nuclear properties, the calculation of neutrino-nucleus cross 
sections is essentially parameter free. In the particular model introduced in this work, once 
the ground state of a target nucleus is calculated with the RHB model, the only additional 
parameter in the PN-QRPA calculations of excited states is the strength g' of the zero-range 
Landau-Migdal term, which accounts for the contact part of the nucleon-nucleon interaction. 
The value of this parameter depends on the isovector properties (symmetry energy) of the 
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effective interaction in tlie particle-fiole cliannel, and is adjusted to reproduce tlie excitation 
energy of the GTR in ^°^Pb. The same value of g' is then used in calculations of charge- 
exchange excitations in all other nuclei. Since it does not contain parameters adjusted to 
specific target nuclei, or even to a particular mass region, the model can be extended to 
regions of nuclei far from stability for which no data on ground and/or excited states are 
available, including those on the path of the r-process. 

We have performed a number of illustrative test calculations for charged-current neutrino 
reactions on ^^C, ^^O, ^^Fe, and ^°^Pb, in the low-energy range below 100 MeV neutrino en- 
ergy The results have been compared with those obtained in previous theoretical studies 
based on the shell model and the non-relativistic random phase approximation, and with 
the available experimental values for flux-averaged cross sections. In addition to the total 
neutrino-nucleus cross sections, we have also analyzed the evolution of the contributions 
of different multipole excitations as a function of neutrino energy. It has been shown that 
except at relatively low neutrino energies < 30 MeV, for which the reactions are domi- 
nated by transitions to IAS and GTR states, at higher energies the inclusion of spin-dipole 
transitions, and also excitations of higher multipolarities, is essential for a quantitative de- 
scription of neutrino-nucleus cross sections. Finally, we have also investigated the cross 
sections for reactions of supernova neutrinos on ^^O and ^ospi-, target nuclei as functions of 
the temperature and chemical potential. 

The results for the test cases are in good agreement with the few available data, and with 
the cross sections calculated in the shell model for reactions on hght nuclei. The advantage of 
the RHB plus PN-RQRPA model over the shell model approach is, of course, the possibility 
of performing calculations for reactions on heavier nuclei and in regions of nuclei far from 
stability. The differences between the results of various RPA-based calculations, especially in 
heavier nuclei and at low neutrino energies, can most probably be attributed to the different 
effective interactions used in modeling the structure of target nuclei, but they also indicate 
that more detailed experimental and theoretical studies of the transitions that contribute to 
the neutrino reaction rates must be performed. 

After the present test study, future applications of the RHB plus PN-RQRPA model will 
include a more extensive investigation of charged-current neutrino-nucleus reaction rates, as 
well as processes induced by neutral-current neutrino reactions, and neutrino nucleosynthe- 
sis. 
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FIG. 1: Inclusive neutrino-nucleus cross sections for the ^'^'^Pb(z^e! e^)^'^^Bi reaction, evaluated 
using the effective momentum approximation (dashed) , the Fermi function correction (dotted) , and 
their combination (solid). The PN-RQRPA calculation is performed with the DD-ME2 effective 
interaction. 
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FIG. 2: The RHB plus PN-RQRPA inclusive neutrino-nucleus cross sections for the ^^C{V(,, e^)^^N 
reaction. The different curves correspond to cross sections evaluated by successively increasing the 
maximal allowed angular momentum in the sum over J in Eq. from Jmax = 0^ to Jmax = 7^. 
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FIG. 3: Inclusive cross sections for the i^O(z^e, e~)i^F, 4"Ar(i/e, e-)4"K, and ^"^Pb(z^e,e 



\208 



Bi 



reactions. The RHB plus PN-RQRPA(DD-ME2) results are compared with those obtained in the 



QRPA with the Skyrme interaction SIII 
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FIG. 4: Inclusive cross section for the ^^Fe(t'e, e^)^^Co reaction. The solid curve corresponds to 
the full calculation including all transitions with multipolarity J < 7^, whereas the cross section 
evaluated for the IAS and GT transitions only is represented by the dashed curve. The dotted 
curve is the difference between the two calculations. 
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FIG. 5: Contributions of different multipole transitions to the RHB plus PN-RQRPA inclusive 
cross section for the ^^Fe(z/e, e~)^^Co reaction, at the incoming electron neutrino energies E,y^ = 
20, 40, 60, and 80 MeV. 
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FIG. 6: Inclusive cross section for the O(z^e) e ) F reaction, averaged over the supernova neutrino 



flux and plotted as a function of temperature. The PN-RRPA results obtained wit 



h the DD-ME2 



eff'ective interaction are compared with the shell-model cross-sections from Ref. 
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FIG. 8: Contour plot of supernova neutrino-nucleus cross sections for the ^^0(fe, e~)^^F reaction, 
as functions of the temperature T and chemical potential a which determine the Fermi-Dirac 
neutrino spectrum Eq. (j38p . The solid curves in the plot correspond to constant values of the cross 
section, and lighter shading denotes the increase of cross sections. 
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TABLE I: Flux-averaged neutrino- nucleus cross sections for the C target nucleus. The results of 
several shell- model and QRPA based calculations, including the present PN-RQRPA(DD-ME2), 
are compared with available data. 







{^e,e ) 




(cj) (10-^0 cm2) 


(a) (10-42 


RPA(SIII) [22] 


19.23 


55.10 


QRPA(SIII) [22j 


20.29 


52.0 


CRPA(WS+LM) [29] 


18.18 


19.28 


PQRPA(PIII) [38] 




17.54 


RPA [58] 




13.60 


RPA(SGII) [20] 


13.5 


12.9 


RPA(SIII) [20] 


14.5 


16.5 


SM(WS(0+l+2+3)^w) 


[19] 13.2 


12.3 


SM(On(j xO.64) [56] 


19.2 


15.1 


SM(HF(0+l+2+3);ia;) 


[22] 15.18 


16.42 


PN-RQRPA(DD-ME2) 


19.59 


12.14 


Exp.(KARMEN) [8] 




14.0±1.2 


Exp.(LSND) [54, 55] 


12.4±0.3±1.8 


13.2±0.4±0.6 


Exp.(LAMPF) [10] 




14.1±2.3 
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TABLE II: Flux-averaged cross sections for the v^, reaction on O, Fe, and Pb target nuclei. 



i*^0(i/e, e~)i'^F 56Fe(zy„ e-f^Qo 2°^Pb(z^e, e^ f^^^l 
(ct) (10-42 gm^) (a) (10-42 gm^) (a) (10-^2 cm^) 



SM(Ona; xO.64) [56] 10.8 

TM [591 214 



TDA(SKIII) [60] 2954,3204 
RPA [581 14.55 277 2643 



CRPA(WS+LM) [6]J 240 3620 

(Q)RPA(SIII,SGII) 16.90,17.20 [20] 352 [57] 4439 [23 

PN-RQRPA(DD-ME2) 13.18 140 2812 

Exp.(KARMEN) [9, 61] 256±108±43 
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